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We study the enumeration of alternating links and tangles, considered up to topo- 
logical (fiype) equivalences. A weight n is given to each connected component, and in 
particular the limit n yields information about (alternating) knots. Using a finite 
renormalization scheme for an associated matrix model, we first reduce the task to that of 
enumerating planar tetravalent diagrams with two types of vertices (self-intersections and 
tangencies), where now the subtle issue of topological equivalences has been eliminated. 
The number of such diagrams with p vertices scales as 12^ for p ^ oo. We next show 
how to efficiently enumerate these diagrams (in time ~ 2.7^) by using a transfer matrix 
method. We give results for various generating functions up to 22 crossings. We then 
comment on their large-order asymptotic behavior. 
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1. Introduction 

The mathematical description of knots and related combinatorial objects is a subject 
with a history of several centuries. Giving an account of the progress made over the 
years is beyond the scope of the present paper, and we refer to [Q for an introduction to 
the existing literature (from the 19th century onwards). An old idea, first published by 
Ozanam in 1694, is to treat a knot as a planar diagram with over- and underpasses. 
This trick, however simple it may appear, is essential to the present work, since it allows 
us to effectively reduce the dimension of the problem from three to two, and in particular 
to apply the powerful tools of matrix models in the planar limit. The classification of 
such diagrams with respect to topological equivalences calls for a particular geometrical 
framework (tagged geometria situs by Leibniz 0), as was probably first pointed out by 
Vandermonde Q in 1771. 

For technical reasons, we shall here constrain our attention to the subset of alternating 
knots (and various generalizations such as links and tangles), which by definition possess 
a representation as a diagram in which over- and underpasses alternate. The first efforts 
in this domain date from the late 19th century through the pioneering works of Tait, 
Kirkman and Little. Tait originally thought that all knots were alternating, but Little 
made it evident (though no proof was known at that time) that non-alternating knots 
start appearing at order p = 8 intersections. It is now believed that alternating knots 
are asymptotically subdominant (the corresponding statement has been proven rigorously 
in the case of links 0). However, the enumeration of alternating knots remains a very 
appealing mathematical problem, and to this date no explicit formula giving the number 
of knots at order p is known. 

The subtle question of factoring out topological equivalences for alternating knots 
is facilitated by two remarkable theorems which were first conjectured by Tait, but only 
proved rigorously quite recently. First, alternating diagrams can be reduced to diagrams 
having minimal crossing number by iteratively eliminating irrelevant crossings through the 
move shown in Fig. 3 a |^ . Second, two such reduced alternating diagrams are topologically 
equivalent if and only if they are connected by a sequence of simple transformations (see 
Fig. 4) known as flypes . 

In the first part of the paper we shall show that counting alternating knots modulo 
the flype equivalence is equivalent to computing the generating function for a certain class 
of planar tetravalent diagrams that are illustrated in Fig. 2. These diagrams possess two 
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types of vertices (self-intersections and tangencies), which are weighted separately. The 
proof is inspired by the renormalization procedure originally invented within the context of 
quantum field theory (QFT), and here applied to the random matrix model that generates 
the above-mentioned planar diagrams. In particular we show that factoring out topological 
equivalences is tantamount to including various counterterms in the effective action of the 
matrix model. Although we have found it appropriate to state the argument using the 
language of physics, it should be stressed that it eventually involves nothing but formal 
manipulations of generating functions, and thus is completely rigorous. 

Unfortunately we have not found a way of analytically computing the generating 
function of the diagrams depicted in Fig. 2 for general n, although a number of cases 
{n = 1 [||,n = 2 0,n = cxD and n = —2 [jl^) can be solved exactly. In the second part 
of the paper we therefore proceed to describe how the integer coefficients of the generating 
function can be found, up to any desired order p, using a computer. The number of 
diagrams to be enumerated scales as 12^ for p ^ 1, but our algorithm does the counting 
in a most expeditious fashion, only using a time ~ 2.7^. Once again, we have inspired 
ourselves from theoretical physics, reformulating the counting process as the action of a 
transfer matrix, which is an object normally used to describe the discrete time evolution of 
a quantum system. Mathematicians may think of the transfer matrix as a linear operator 
that acts on a set of basis states consisting of appropriately defined sub-diagrams. 

From the description given below it will be evident that the applicability of our transfer 
matrix method extends far beyond the context of the present paper. Indeed, it can be 
generalized to the computation of generating functions for all sorts of planar diagrams, 
provided that the latter have at least one external leg. In particular, we are presently 
unable to generate results for closed knots (usually known in the mathematics literature 
simply as knots, as opposed to tangles which possess external legs) . 

We finally state our results, here given up to a maximal order ofp = 22, and we analyze 
the asymptotic behavior of the various generating functions. The number of objects at 
order p ^ 1 is generally supposed to have the asymptotic behavior C/U^p~", where, in 
the language of statistical physics, /U is a critical temperature and a a critical exponent. 
Amazingly, the theory of two-dimensional quantum gravity provides a conjecture for a, 
based solely on the 0(7i) symmetry of the underlying matrix model, whereas the leading 
scaling exponent fi can only be accounted for in the exactly solvable cases. Comparing 
the conjecture for ct(n) with our numerical results, we find good agreement in its expected 
range of validity (0 < n < 2). 
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2. The matrix model and its renormalization 

The objects we want to consider are tangles with 2k "external legs", that is roughly 
speaking the data of k intervals embedded in a ball B and whose endpoints are given 
distinct points on the boundary dB, plus an arbitrary number of (unoriented) circles em- 
bedded in B, all intertwined, and considered up to orientation preserving homeomorphisms 
of B that reduce to the identity on dB. Tangles with 4 external legs will be simply called 
tangles. As mentioned in the introduction, we represent these objects using diagrams, that 
is regular projections on a plane, and restrict ourselves to alternating diagrams for which 
under- and over-passes alternate as one follows any connected component. 

A natural framework is to introduce a matrix model which generates in its Feynman 
diagram expansion these link diagrams. We shall define it now. 



2.1. Definition of the 0{n) matrix model 



As in ||Tl|,p|,pI]|| , we start with the following matrix integral over N x N hermitean 
matrices 



(2.1) 



a=l 



where n is a positive integer. 

Expanding in power series in g generates Feynman diagrams with double edges ( "fat 
graphs") drawn in n colors in such a way that colors cross each other at the vertices. 
By taking the large N limit one selects the planar diagrams, which can be redrawn as 
alternating link diagrams, cf Fig. 1. 





Fig. 1: A planar Feynman diagram of (^T|) and the corresponding alternat- 
ing link diagram. 
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More precisely, the large "free energy" 

logZW(n,ff) 

^^^'^) = il^oo m ^^-^^ 

is a double generating function of the number fk-p of alternating link diagrams with k 
connected components and p crossings (weighted by the inverse of their symmetry factor, 
and with mirror images identified): 

oo oo 

^(^'^?) = EE/fc;^'^V (2.3) 

k=l p=l 

If one is interested in counting objects with a weight of 1, one cannot consider the 
free energy which corresponds to closed diagrams, but instead correlation functions of the 
model which generate diagrams with external legs, that is tangle diagrams. Typically, we 
shall be interested in the two-point function 

G(n,,) = ^lim^(ltrM„2^ (2.4) 

where the measure on the matrices Ma is given by Eq. (|2.1|) and a is any fixed index, which 
generates tangle diagrams with two external legs; and the connected four-point functions 

ri(n,(7)= lim /i-tr(M„M,)A (2.5.1) 

r^{n,g) = hm [\-MMlMl^-G(n,gf (2.5.2) 



where a and h are two distinct indices, which generate tangle diagrams with four external 
legs of type 1 and 2 (see Fig. 2). 
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2.2. Renormalization of the 0{n) model 

The model presented above counts tangle diagrams, but not tangles. There is a re- 
dundancy in the counting since to a given tangle/link corresponds many diagrams. In 
order to properly count tangles, this redundancy must be removed. In the case of alter- 
nating diagrams one can distinguish two steps. First one must find a way to select only 
reduced diagrams which contain no irrelevant crossings (Fig. 3a)); such diagrams will have 
minimum number of crossings. It turns out to be convenient to introduce at this point 
a closely related notion: a link is said to be prime if it cannot be decomposed into two 
pieces in the way depicted on Fig. 3 b). It is clear that at the level of diagrams, forbidding 
decompositions of the type of Fig. 3 b) automatically implies that the diagram is reduced; 
and we shall therefore restrict ourselves to prime links and tangles. 





b) 

Fig. 3: a) An irrelevant crossing, b) A non- prime link. 

There may still be several reduced diagrams corresponding to the same link: according 
to the fiyping conjecture, proved in , two such diagrams are related by a finite sequence 
of fiypes, see Fig. 4. 





Fig. 4: A fiype. 

To summarize, there are two problems: a) the diagrams generated by applying Feyn- 
man rules are not necessarily reduced or prime; b) several reduced diagrams may cor- 
respond to the same knot due to the fiyping equivalence. Following the study of |1T0[] , 



we note (cf Figs. 3 and 4) that this "overcounting" is local in the diagrams in the sense 
that problem a) is related to the existence of sub-diagrams with 2 external legs, whereas 
problem b) is related to a certain class of sub-diagrams with 4 external legs. Clearly such 
graphs can be cancelled by the inclusion of appropriate counterterms in the action. We 
are therefore led to the conclusion that we must renormalize the quadratic and quartic 
interactions of Renormalization theory tell us that we should include in the action 



from the start every term compatible with the symmetries of the model, since they will be 
generated dynamically by the renormalization. A key observation is that, while there is 
only one quadratic 0(?T.)-invariant term, there are two quartic 0(n)-invariant terms, which 
leads to a generalized model with 3 coupling constants in the action: 

Z"- '{n,t,gi,g2)= [[dMae V 

a=l 

(2.6) 

The Feynman rules of this model now allow loops of different colors to "avoid" each other, 
which one can imagine as tangencies (Fig. 5). 



Si 



92 



Fig. 5: Vertices of the generalized 0{n) matrix model. 



We define again the correlation functions G{n,t, gi, §2) and ri{n,t, gi, §2) (Eqs. ( PT^ ) 
and (2.5)), and want to extract from them the counting of colored alternating tangles with 
external legs. 

The idea is to find the expressions of t{g), gi{g) and g2{g) as a function of the renor- 
malized coupling constant g, in such a way that the overcounting is suppressed and the 
correlation functions are generating series in g (and n) of the number of colored tangles. 



The derivation of the renormalization equations was performed in [|T0[ . Consider first 
the removal of irrelevant crossings and non-prime links. It is clear that one must remove 
all two-legged subdiagrams, that is impose 



G{n,t,gi,g2) = 1 



(2.7) 



This implicitly fixes t{g). 






Fig. 6: Breaking a flype into two elementary flypes. 
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Next, consider the flyping equivalence. A flype can be made of several "elementary" 
flypes (Fig. 6), an elementary flype being by definition one that cannot be decomposed any 
more in this way. In the terminology of QFT, these elementary flypes are one simple vertex 
connected by two edges to a non-trivial two-particle irreducible in the horizontal channel 
tangle diagram (H-2PI). We need to introduce auxiliary generating functions H[{g), -^2(6') 
and V2{g) for non-trivial H-2PI tangles of type 1, of type 2 and of type 2 rotated by 7r/2 
respectively. By considering all possible insertions of elementary flypes as tangle sub- 
diagrams of a diagram, we flnd (Fig. 7) that the renormalization of gi and g2 is given 
by: 

gi{g) = gil~2H!,{g)) (2.8.1) 
92{9) = -9{H[{g) + V^{g)) (2.8.2) 
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Fig. 7 : Counterterms needed to cancel flypes. 



All that is left is to flnd the expressions of the auxiliary generating functions in terms 
of known quantities. They are easily obtained by decomposing the four-point functions in 
the horizontal and vertical channels 

1 



H'±H[ = 1- — 
2 1 ^ 



1 

(i-^)(i + (n-M)r2 + ri) 



(2.9a) 
(2.9b) 



The three renormahzation equations (|2.7| ) and (2.8), supplemented by the expressions (2.9), 
are enough to fix t{g), gi{g) and 02(9) ■ Solving them gives access to the F^, which are the 
generating series of the numbers of prime alternating tangles of type i. However, we can 
go further. By computing other correlation functions in the model and composing them 
with the solutions t{g), gi{g)^ 92(9) of the equations above, one can extract the generating 
functions of the number of alternating tangles with an arbitrary number of external legs. 
An example will be given below. 



3. An algorithm for the counting of decorated planar diagrams 

Unfortunately the matrix model described in section 2.1 cannot be solved for a generic 



value of n (see however [p!2|,p|,p!0[| for solutions of particular values). We therefore turn to 
the description of an algorithm that allows to enumerate the required diagrams to a given 
order p of crossings. This is one of several similar algorithms described in [1T3|,|9[|. The 



numbers of diagrams thus computed are the input needed for the combinatorial treatment 
of section 2.2. 

The principle of the algorithm is to iterate a transfer matrix which "builds" the 
diagrams slice by slice: starting from an initial state consisting of all external legs (two 
in the case under consideration), the system is time evolved through the addition of p 
intersections, until an empty final state is obtained. 

We shall first concentrate on the enumeration of (prime, alternating) tangle diagrams, 
which correspond to the "unrenormalized" model of section 2. Adding tangencies, which is 
needed to take into account the fiyping equivalence, will be discussed in Section 3.2, since 
it is an elementary extension of the algorithm. 

3.1. The single-step algorithm 

The algorithm we want to describe, called "single-step algorithm" in 0, is best ex- 
plained by an example. Let us consider the tangle diagram shown in Fig. 8. Let us start 
from an initial state given by the two external legs (edges A and G) . We want to "evolve" 
the various open lines by following them and adding the new crossings they meet. 

Moving along either of the edges A or G, a new line segment (DE resp. HI) is encoun- 
tered. The question then arises which of these to process first. We resolve this ambiguity 
by stipulating that in any given state, we evolve the line which at that instant is uppermost. 
At time t = 1, the edge A thus becomes B, and the new line segment DE is added. The 
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Fig. 8: Working principle of the single-step algorithm, a) A two-legged 
knot diagram with p = 6 intersections and = 3 connected components. 
The edges are labelled from A to M. b) The same diagram in the time- 
slice representation. For reasons of clarity, the time slices are not drawn in 
chronological order. 

edge D is now the new top line. Analogously, at the instants t = 2 and t = 3, the top line 
(D resp. K) crosses a new line segment, which is then added to the current state. We can 
formalize this by stating the transformation rule shown in Fig. 9.1. 

At time t = 3, the new top line carries the label B, which was however already 
produced by the transformation acting at t = 1. We therefore proceed, at t = 4, to the 
identification of the two "copies" of B, joining them through an arch. This is an example 
of the general transformation rule shown in Fig. 9.2. The addition of an arch means that 
the lines intermediate between the two instances of B (at positions p = 1 and q on Fig. 9.2) 
can henceforth not communicate with the lines at the exterior of the arch. These "trapped 
lines" must therefore eventually evolve to the empty state (vacuum), independently of the 
rest of the diagram. Since the current top line must always be treated first, the evolution of 
a possible set of "trapped lines" must take place at a later time. Thus, on Fig. 8, we cannot 
always draw the time-slices in chronological order. The existence of a number of trapped 
lines is visualized on Fig. 9.2 by a delimiter (shown as a gray rectangle), which separated 
the remaining lines into two blocks. Lines in different blocks cannot communicate, and so 
the transformation 2) only applies when p and q belong to the same block. 
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Fig. 9: The two types of transformations in the single-step algorithm. 1) 
Addition of a new line segment. 2) Identification of the top line (at position 
p) with another line (at position g), accompanied with the creation of a new 
block. 

The purpose of the transfer matrix is not only to count the total number of tangle 
diagrams, but to do so for any fixed number of connected components. In particular, when 
performing a type 2 transformation, we need to know whether the points p and q were 
already connected though an arbitrary number of edges at an earlier time. On Fig. 9 we 
have represented this information by a number of lines on the left, connecting the points 
at a given instant into pairs. It may thus happen that on Fig. 9.2, r = p and s = q. In 
this case, the type 2 transformation marks the completion of one connected component in 
the tangle diagram. 



10 11 12 13 



Fig. 10: Intermediate states produced by applying the single-step algorithm 
to the tangle diagram shown in Fig. 8. 



We are now ready to define the set of states on which the transfer matrix acts. A 
state is defined by an even number of points (represented on Fig. 9 as black rectangles), 
connected into pairs by means of the edges encountered at previous times. In addition, 
the points are divided into £+1 blocks by means of £ > delimiters. On Fig. 10 we show 
the set of intermediate states corresponding to the time-slice representation of the tangle 
of Fig. 8. 
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Finally, we must define the transfer matrix T which counts all tangle diagrams with p 
vertices and k connected components. Its entries T^^, where a and b are two basis states of 
the kind just defined, are unless 6 is a descendant of a. An allowed state 6 is a descendant 
of a if it can be obtained via a transformation of one of the two types shown on Fig. 9 (for 
an arbitrary even q > 2 belonging to the same block as p = 1), followed by an arbitrary 
number of simplifications. Tab is then the sum over all transformations from a to 6 of the 
corresponding weight: 1 or n depending on whether one closes a connected component or 
not. The simplications are elimination of superfiuous delimiters, and reduction of states 
to a "normal form" by using various symmetries of the transfer matrix [^. 

3.2. Tangencies 

Until now we have been discussing the enumeration of tangle diagrams in which every 
vertex represents a crossing. However, to account for the flype equivalence we need to 
enumerate more general diagrams with pi intersections and p2 tangencies, as discussed in 
Section 2. Adding a tangency rather than an intersection is obtained by modifying the 
transformation if Fig. 9.1, so that the two points added at time t + 1 are both immediately 
above (resp. immediately below) the uppermost point at time t. any given time, specifies 
how many tangency transformations were used prior to that instant. The desired diagrams 
are then generated by sequences of pi crossing transformations (type 1), p2 trangency 
transformations, and pi + p2 + 1 transformations of type 2, so that no intermediate state 
is empty. 



4. Numerical results and analysis 

The implementation of the above algorithm provides an efficient numerical way to 
count the diagrams needed for our purposes. We present here some of the data obtained 

mi- 



m 



First, we give one example of the "raw data" obtained from the algorithm. It cor- 
responds to counting a certain set of diagrams without taking into account topological 
equivalences: in the present case, connected two-legged alternating diagrams (called two- 



legged knot diagrams in |]T^), see Tab. 1. The 2PI diagrams are precisely the reduced 
diagrams of prime tangles. The importance of going to orders as high as possible comes 
from the issue of the asymptotic behavior which will be discussed below. 
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Tab. 1: Total number G of two-legged diagrams, as well as the subset E2 of 
2PI diagrams. 

Let us also give two examples of renormalized data. Tab. 2 shows the number of four- 
legged prime alternating tangles of the two possible types (Fig. 2) with arbitrary number 
of connected components, while Tab. 3 shows the number of six-legged tangles, which exist 
in five types (Fig. 11). 

Finally, let us discuss the expected asymptotic behavior of these series. Let us consider 
an unrenormalized generating series for tangle diagrams, say G{n,g) — J2^i^p{''^)9^- 
Then it is expected that 

ap(n)^^°°e^^"^^j9-"(") (4.1) 

The conjecture holds for any correlation function in the model; for the free energy F(n, g) = 
Z^^i fp{''^)g^i that is the generating series for link diagrams (the symmetry factors being 
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Tab. 2: Table of the number of prime alternating tangles with 4 external 
legs. 
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Tab. 3: Table of the number of prime alternating tangles with 6 external 
legs. 

most likely irrelevant for the asymptotic behavior), it is 

fp{nf^°°e'^''^Pp-''^''^-^ (4.2) 
13 








Fig. 11: The five types of tangles with 6 external legs. 

The bulk entropy s(n) can be easily extracted with a good accuracy from the numerical 
data. However it is non universal, which implies that it is not preserved by the renormal- 
ization of the model. On the other hand, the exponent a{n) is universal, and presumably 
preserved by the renormalization (unless an unusual phase transition takes place in the 
model). Since the raw data are "cleaner" than the renormalized data, they are better 
suited to the extraction of a(7i). A conjecture made in [rT|, based in particular on the 
KPZ relation [14|, is that for n analytically continued to < n < 2, one has 



a(n = -2cos(7rzy)) = 1 + 1/zy l/2<p<l 



(4.3) 



Using Tab. 1 one can try to evaluate a = a{n = 0), the exponent of tangles with 
minimum number of components, which also gives the exponent of knots according to 
Eq. (^]2|). The conjecture ( ^?3|) yields a = 3. Unfortunately, with the current data, it 
is difficult to estimate a without any knowledge of the subleading corrections. Indeed, a 
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direct fit gives a ~ 2.76 but with very weak accuracy. In |13], it was suggested to fit the 
data with = e*Pp~°(alogp + 6 + the result being in good agreement with the 

conjecture: 

e' = 11.416 ±0.005 a = 2.97 ±0.06 a = 0.04 ± 0.02 6 = 0.1 ±0.03 (4.4) 

However there is httle evidence for such a logarithmic correction, and the issue remains 
open. A more detailed discussion of a{n) for all n can be found in ||]. 
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